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ABSTRACT 

We report the detection of the 62 /xm feature of crystalline water ice in emission towards the bow- 
shaped Herbig-Haro object HH 7. Significant amounts of far infrared continuum emission are also de- 
tected between 10 and 200 ^m, so that Herbig-Haro objects cease to be pure emission-Une objects at 
FIR wavelengths. The formation of crystalline water ice mantles requires grain temperatures Tg^ ^ 100 
K at the time of mantle formation, suggesting that we are seeing material processed by the HH 7 shock 
front. The deduced ice mass is ~ 2 x 10~^ M© corresponding to a water column density N(H2 0)~ 10^^ 



cm^ ; an estimate of the [H20]/[H] abundance yields values close to the interstellar gas-phase oxygen 
abundance. The relatively high dust temperature and the copious amounts of gas-phase water needed 
to produce the observed quantity of crystalline water ice, suggest a scenario where both dissociative and 
non-dissociative shocks co-exist. The timescale for ice mantle formation is of the order of ^ 400 years, 
so that the importance of gas-phase water cooling as a shock diagnostic may be greatly diminished. 

Subject headings: (ISM:) dust, extinction — ISM: Herbig-Haro objects — ISM: individual (HH7) — 
ISM: lines and bands — infrared: ISM: continuum — infrared: ISM: lines and bands 



1. INTRODUCTION 



Herbig-Haro objects (HH; Haro 195C , Herbig 1951 ) are 
emission-line objects acting as signp osts fo r the shock re- 
gions (Draine, Roberge & Dalgarno 1983 ; HoUenbach & 
McKee |l989 ) originating at the interface between stellar 
winds accelerated by Young Stellar Objects (YSOs) and 
the circumstellar or cloud ambient material. The temper- 
ature of dust grains in these shock regions can raise to only 
a few hundred degr ees at most in dissociative shocks (Hol- 
lenbach & McKee 198£ ), so that their thermal emission 
cannot be detected below 10 /xm. 

The instruments on board the Infrar ed Space Observa- 
tory satellite (ISO, Kessler et al. 1996| ) opened unprece- 
dented possibilities for far infrared continuum studies of 
cold objects, including HH objects. In this Letter we 
presen t the data obtained with the Long (LWS, Clegg 



trum for HH 7 are discussed. Nomenclature for the ten 
LWS detectors is described in the ISO Data User Man- 



ual ( http://www. iso.vilspa. esq. es/manuals/lws^idumt ; de- 
tcctors SWl to SW5 (43-90 ^m), and LWl to LW5 (80- 
197 /im), are sometimes referred to as "short" and "long" 
wavelength detectors respectively. 

2. RESULTS 

The LWS beam centered on HH 7, also includes object 
HH 8 somewhat 20" off-axis and HH 10 at the edge of the 
beam; these object s are however fainter than HH 7 at 2 
/xm (Molinari et al. 1999) and the beam profile suppresses 
their possible contribution even more. Apart from a con- 
tamination by the strong near by s ource SVS 13, which will 
be discussed in detail in Sect. 2.2, it is plausible to assume 



et al. |1996|) and Short (SWS, de Graauw et al. |l996D 



Wavelength Spectrometer towards HH 7, the leading bow- 
shaped shock of the HH7-11 chain emanating from the 
YSO SVS 13, in the star forming region NGC 1333 in 
Perseus (d=350 pc). Details about the observations 
an d dat a reduction are given elsewhere (Molinari et 



that HH 7 dominates the observed spectrum. 
2.1. The 62 ^m feature 



In Fig. 2.1 



al. 1999| ). In Sect, g the additional data analysis pro- 



cedures we adopted to derive a reliable continuum spec- 

"^Based on observations with ISO, an ESA project with instruments 
Germany, the Netherlands and the United Kingdom) with the participat 



the complete spectra observed towards HH 7 
and SVS 13 are shown. A broad feature extending from 
roughly 50 to 70 /J,m is clearly visible in the HH7 spec- 
trum. For clarity the SW3 detector spectrum, with an 
offset applied to align it with the continuum of the ad- 
jacent detectors, is shown with plus symbols. The pri- 
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mary concern was to make sure that the 50-70/iin fea- 
ture was not a result of a residual instrumental effect and 
that it is intrinsic to HH 7. The passband calibration is 
known to be somewhat inaccurate for detector SWl (it 
may be worse than 50% both in absolute and relative 
terms) which is the most sensitive to transients. Detec- 
tors SW2, SW3 and SW4 however, are stable in relative 
(passband) terms and a conservative figure of 30% can be 
assumed for their absolute calibration accur acy . The possi- 
bility that the broad feature visible on Fig. 2.1 may be due 
to the near-IR leaks present in the LWS detecto rs filters 



{http://isowww.estec.esa.nl/notes/lws_0197.htm.) can also 



be excluded. Finally, this feature might result from con- 
tamination effects from the bright nearby source SVS 13, 
which is the candidate exciting source for HH 7; these ef- 



fects will be discussed in detail in the Sect. 2.2 
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Fig. 1 - LWS averaged spectra observed towards HH 7 (bot- 
tom) and SVS 13 (top). The different portions are the 10 LWS 
detectors, labeled with the names used in the text. The plus 
symbols for HH 7 represent detector SW3 with an applied offset 
of 5 Jy to bring it in line with the adjacent detectors. 

Here we just note here that if the 50-70 /xm feature 
seen on the HH 7 spectrum was due to a fraction of the 
flux emitted by the nearby contaminating source SVS 13, 
we would expect to see the same feature at a comparable 
"line-to-continuum" ratio on the SVS 13 continuum spec- 
trum; Fig. 2.1 clearly shows that this is not the case. We 
conclude that the 50-70 /.tm feature is real and intrinsic to 
HH 7 and we consider it as a 2cr detection, assuming that 
the noise of the feature is equal to half the gap between 
detectors SW3 and SW2. We identify this as the 62/xm fea- 
ture due to the longitudinal acoustic mode s of crystalline 
water ice (Bertie, Labbe & Wha lley | l969| ), observed for 
the first time by Omont et al. ( 19901) in the expanding 
envelopes of post-AGB stars; ISO-LWS has det ected this 
feature in the spectra of similar objects (Barlow 1998) and 



of a few Herb ig Ae /Be stars (Waters & Waelkens 



1998 



Malfait et al. 199S| ). It is the first time this feature has 



been detected towards Herbig-Haro objects. 

2.2. The 2-200 fim Continuum 

Several SWS line scans were used to estimate the contin- 
uum at different wavelengths for A <40 /im. The line-free 
portions of each individual spectrum were used to build 



flux histograms, and a gaussian was then fitted to the core 
of the distribution to obtain the centroid (average fiux) 
and the standard deviation; the latter was then divided 
by the square root of the number of points in histogram 
to get an estimate of the uncertainty. These uncertainties 
reflect the internal accuracy of the estimates; infact the 
true uncertainty may be higher. In particular, all data 
shortward of ~ 10 /im are at the detection limit of the 
SWS and they will be treated as Icr upper limits. 

LWS scans were averaged using the ISO Spectral Analy- 

sis Package (ISAP, http://www.ipac.caltech.edu/iso/isap/isap.htmi) 
using a median clipping algorithm which rejected outly- 
ing data points due to incomplete removal of glitches 
from cosmic rays hits. However, contamination effects 
in the HH7 spectrum due to the nearby {^ 70" NW) 
candidate exciting source SVS 13 must be taken into ac- 
count to assess if the observed continuum (Fig. 2.1 is 
all due to HH 7; for this purpose we used an irregu- 
larly sampled raster map of Mars taken as part of the 
LWS calibration programme. We first us ed the LIA 
( http://www. ipac. caltech. edu/iso/lia/lia. htm\ } routine IN- 
SPECT_RASTER to determine the relative position of 
HH7 and SVS 13 in the [Y,Z] spacecraft frame of refer- 
ence. For each Mars raster position along the direction 
connecting the two sources, we averaged the Mars spectra 
detector by detector; each detector average at the var- 
ious off-axis positions was then ratioed to the analogue 
detector average of the on-axis spectra, and a set of ten 
contamination factors for each position was obtained. A 
spline interpolation was then used to estimate these fac- 
tors (one per detector) at an off- axis distance of 68", the 
distance between HH7 and SVS 13; we find values 0.024, 
0.025, 0.037, 0.05, 0.075, 0.09, 0.10, 0.11, 0.13 and 0.15, 
which increase with wavelength as expected due to diffrac- 
tion. We multiplied the observed SVS 13 spectra by these 
numbers, and the resultant spectrum was subtracted from 
the observed HH 7 spectrum. 

This method for estimating the contamination from 
nearby sources inevitably suffers from the irregular sam- 
pling of the Mars raster map. Raster points did not lay 
exactly along the HH 7-SVS 13 direction on the focal plane, 
and the interpolation between the correction factors esti- 
mated at each relevant raster position introduced an addi- 
tional uncertainty. At the end of the procedure (see Fig. ||) 
we find worse alignment between adjacent LW detectors, 
and detector SW5 goes to negative flux values (not re- 
ported on Fig. g). We beheve this may be due to the 
fact that as the magnitude of the applied contamination 
correction increases with wavelength, so it does the asso- 
ciated uncertainty; in particular, detector SW5 is in the 
critical region where the contamination correction starts 
to be high (2-3 times higher than for SW2 and SW4) 
while the observed signal is still low (similar to SW2), 
so that a slightly overestimate of the contamination frac- 
tion is enough to bring its corrected values below zero. 
The LWS LW detectors are those for which the contami- 
nation fractions are higher and for which diffraction and 
source's extension effects (which here have been neglected) 
are more severe. An additional complication, which does 
not affect the SW detectors, is that heavy fringing is ob- 
served (which has been here removed using standard tools 
available in ISAP); this has the effect of modulating the 
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beam size as a function of wavelength even within indi- 
vidual detector bands. We did not take this into account, 
instead deriving a single contamination factor per detec- 
tor; multiplying or dividing by a constant, however, has 
the effect of changing the slope. Hence, although the ob- 
served continuum levels are high enough to ensure that an 
important part of the observed signal can confidently be 
assigned to HH7, the exact absolute fluxes and spectral 
shape in this wavelength range remain highly uncertain. 

On the short wavelength side of the LWS spectrum on 
the other hand, the existence of intrinsic continuum emis- 
sion from HH7 depends more critically on the particular 
value of the contamination factors; doubling this factors 
would lower the spectrum to negative flux levels. If, how- 
ever, the A < 80 /Ltm flux observed towards HH 7 were due 
to contamination, we would expect to see the 62 ^m emis- 
sion feature also on the continuum spectrum of SVS 13 and 
with a comparable "line/continuum" ; instead, no trace 
of such a feature is seen on the SVS 13 continuum (see 
Fig. 2.1). We conclude that FIR emission which is intrin- 
sic to HH7 has been detected, showing for the first time 
that Herbig-Haro objects cease to be exclusively emission- 
line objects at FIR wavelengths. 

3. DISCUSSION 

In the following discussion we will assume that the FIR 
and mm continuum arise from the same region of space, 
based on the spatial coincidence between the HH objects 
and the mm emission distribution; curre ntly availab le 50- 
100 /im data on th e region (Harvey et al. 1984, 1998| ; Jen- 
nings et al. 1987) however, do not have enough spatial 



the observed SED; considering the discussion in Sect. 2.2 



resolution or sensitivity to support this claim. 

The spectral energy distribution (SED) from HH 7 will 
be modeled as thermal emission from dust grains com- 
posed of a silicate core and a water ice mantle. Absorp- 
tivities between 2 and 300 /im were computed with Mie 
theory in the formulation of Wickramasinghe ( |1967D , us- 
ing the comp lex re fractive indices for cr ystalli ne water ice 
(Bertie et al. 1969| ) and silicate (Draine 19851 ). Lon gward 
of 300 /^m the silicate absorptivities by Draine ( 1985| ) were 
adopted; inclusion of water ice mantles may steepe n the 
slope of the Qats vs A relationship (Aannestad 1975] ), pos- 
sibly leading to underestimation of the dust mass. Radia- 
tive transfer is approximated with an analytical treatment 
where dust is distributed on a sphere which is characterised 
by radial density and temperature gradients (Noriega- 
Crespo, Garnavich & Molinari 1998). The presence of 
a cold dust clump centered on the location of HH 7 and 
extending over a larger area than the one traced by the 
optical or near-IR emission (Lefioch et al. 1998), justifies 
a treatment where a dust clump is centrally heated by the 
HH 7 shock. Fig. presents the complete SED towards 
HH 7, aft er removalof the contamination from SVS 13 (see 
Sect. O). The LWS detectors SWl to SW4 have been 
rescaled to a common level preserving their original mean 
value; no other rescaling was made to the rest of the LWS 
spectra. The continuum points from the SWS spectra and 
an estimate for the 1.25 mm fiux integrated on the model 
fit area (see below) are also reported; the latter point has 
to be considered an upper limit since it is contaminated 
by SVS 13 just in the same way it happens for the LWS 
observations. We find that it is impossible to globally fit 



we then decided to give higher priority to the short wave- 
length LWS detectors spectra also because these are the 
ones showing the 62 /im feature, and more closely match 
each other (not considering the applied rescaling). 

While the overall appearance of the continuum depends 
on the radius of the dust clump, its mass and tempera- 
ture gradient, the "line/continuum" ratio of the 62 /im 
feature depends on the relative size of the ice mantle with 
respect to the grain core. The 62 /xm water ice feature 
is well fitted by adopting a core radius of 0.07 /zm and 
a total core -I- mantle radius of 0.1 /xm. The model also 
suggests that the short wavelength end of the LWS spec- 
trum can be identified with one side of the 45 /im water ice 
feature (see also Dartois et al. 1998), also in emission; al- 
though detector SWl is the one that experiences the most 
severe problems in calibration-related aspects, the further 
agreement of the model with the A ~38.4 /im SWS point 
makes this a plausible possibility. The importance of the 
detection of the 62 /im feature is that it is specific to crys- 
talline water ice, as opposed to the 45 /tm feature which is 
instead predicted f or bo th crystalline and amorphous ice 
(e.g. Moore et al. |1994 ). Most of the laboratory studies 



(e.g., Smith et al. 1994 and references therein) show that 
ice deposited onto grains a t T < 100 K is in amorphous 



state (but see Moore et al. 1994 who instead obtain crys- 
talline ice in these conditions) al so for deposition rates as 
high as 20 /<m h~^ (Leger et al. 1983| ); the 62/im feature 



is then an important indicator of the thermal history of 
the grains because it appears at T ^ 100 K when an ir- 
reversible transition from amorphous to crystalline state 
commences. 
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Fig. 2 - Complete spectral energy distribution towards HH 7. 
The full lines represent the LWS spectra after applying the con- 
tamination corrections discussed in the text. Detectors SWl 
to SW4 have been rescaled to a common level preserving their 
original mean value. 

The fit plotted in Fig. ra is obtained adopting a dust 
clump radius of 0.06 pc, which approximates the radius of 
the average LWS beam size at a distance of 350 pc, and 
the si ze of the 1.25 mm continuum emission area (Lefloch 
et al. 1998). The density is assumed constant and equal 



to 6x10 '' cm while the temperature varies from '--^lO 
to 200 K with a ~ —0.4 power-law radial gradient. The 
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bolometric luminosity obtained integrating the model fit 
is Ijfit— 3.7 Lq, while integration of the contamination- 
corrected SED (Sect. p,.2\) yields L^ed ^ 4.5 Lq (since 
A < 10 /xm and A= 1.25 mm data are to be considered 
upper limits); this is about a factor 30 higher than the 
cooling via atomic and m olecular lines observed towards 
HH 7 (Molinari et al. |l999D . The fitted model predicts that 
the bulk of the 62 /im feature is emitted by dust at tem- 
peratures T > 30 K concentrated inside a 4"-radius region 
centered on HH 7, a size comparable to that of the optical 
and near-IR emission which traces the shock front. Since 
these grains should have experienced a rise in tempera- 
ture to values <; 100 K in order for ice mantles to be in 
crystalline state, it is plausible that the 62 /im feature 
originates from dust which has been processed by 
the HH 7 shock. The dust mass in the 4"-radius region 
centered on HH7 amounts to 5xl0~^ Mq (the total dust 
mass implied by the model fit is ~ 0.035 Mq); the relative 
proportion of core and mantle (the core has 70% of the 
total grain radius) implies a water ice mass of 2 x 10~^ 
Mq, or a H2O column density ~ lJ_x 10^* cm^^. H2 
pure rotational lines (Molinari et al. 1999| ) suggest N(H2) 
~ 4.4 X 10^" cm~^ in the same 4"-radius region, imply- 
ing [H20]/[H]~ 1.25 X 10~^, or a factor ~ 4 highe r tha n 
the interstellar O gas-phase abundance (Meyer et al. 19981 ). 
However, this number should be regarded as an upper limit 
because our model assumes that all dust grains are coated 
with ice mantles, which is not necessarily true; besides, 
we cannot exclude the presence of cold H2 ( T ^ 100 K) 



which our ISO observations (Molinari et al. 1999) would 
not trace. This water abundance, even if considered only 
as an order-of-magnitude estimate, is however much higher 
than the gas phase water abundance ( [H2 ]/[H] < 10^^) 
deduced from FIR lines (Molinari et al. |l999| ), and it would 
essentially imply that most of the oxygen is locked into wa- 
ter ice. 

Due to the uncertainty about the 45 /xm feature (see 
above), it hard to tell from the observational viewpoint 
whether gas-phase water produced behind the HH 7 shock 
front (Kaufman & Neufeld 1996| ) is deposited onto bare 
warm grains, or pre-existing ice mantles are warmed-up 
during the passage of a relatively gentle shock f ront. The 
ice optical constants that we used (Bertie et al. 1969 ) are 
from laboratory samples obtained by direct deposition at 
T=173 K and subsequent cooling to 100 K; hence the 
good simultaneous fit of the 62 and 45/im features would 
tend to support the first scenario. The physical condition 
behind low velocity (v^ ^ 40 km s^^), non-dissociative, 
shocks are favourable for the rapid gas-phase incorpora- 



tion o f atom ic oxygen into water (Drai ne, Roberge & Dal- 
garno 1983| , Kaufman & Neufeld 1996). Su bsequent fre ez- 
ing onto grains (Bergin, Neufeld & Melnick 1998, 1999| ) in 
the cooling post-shock region could produce the observed 
water ice mantles, also explaining the minor role played by 
gas-p hase w ater in the cooling of the HH 7 shock (Molinari 



et al. 1999 ). However, non-dissociative shocks are unable 
to raise the grain te mpera ture to the T ^ 100 K (Draine, 
Roberge & Dalgarno 1983) needed to explain the presence 
of crystalline ice mantles; a dissociative componen t, who se 
generated intense UV field (HoUenbach & McKee 1989) is 
much more efficient in heating the grains, is needed. The 
co-existence of both types of shock is indeed a n expe cted 
feature of bow shocks like HH 7 (Smith & Brand 1990) and 



it is indep enden tly supported by FIR lines studies (Moli- 



nari et al. 1999). However, we point out that this relies 



on the general result that crystalline ice requires relatively 
high temperatures to form; if crystalline ice can al so form 
at low (T< 20 K) temperature (Moore et al. 1994) then a 
dissociative shock component may not be needed. 

The alternative possibility that pre-existing amorphous 
water ice mantles are heated up during the passage of the 
shock front seems unlikely since the mantles are easily de- 
stroyed by grain-grain collisions once the shock v elociti es 
exceed ~ 15 km s^^ (Caselli, Hartquist & Havnes 1997). 

We can derive an order-of-magnitude estimate for the 
timescale r^ce of the formation of water ice mantles, divid- 
ing the linear size of the post-shock region by the shock 
velocity; assuming a maximum shock ve locity of '^ 40 



km s 
et al. 



^ for the HH 7 shock (Solf 
1999) and a linear extent 



Bohm 
of - 10' 



1987], Molinari 
the 



for 



HH7 



post-shock regi on as estimat ed fro m H2 2.12/im i mage s 



(Garden et al. [199C 
we estimate Tice 



Everett |1997| , Mohnari et al |1999| ), 
400 yrs, a factor ^ 250 less compared 



to theoretical predictions for ice mantles formation behind 
non-dissociative shocks (Bergin, Neufeld & Melnick 1999| ). 
This very short timescale for ice mantle formation would 
make the gas-phase water cooling practically irrelevant as 
a shock diagnostic. 
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